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METHOD FOR ON-DEMAND GENERATION 
OF INDIVIDUAL RANDOM NUMBERS OF A 
SEQUENCE OF RANDOM NUMBERS OF A 
1/F NOISE 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application is a continuation of copending Intema- 
tioiial Application No PCT/DEOl/04376, Hied Nov. 22, 
2001 , which designated the United States and was cot pub- 
lished in English. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The invention relates to a method for generating sequences 
of random numbers of a 1/f noise. 

Random numbers of a 1/f noise can be used, for example, 
in a transient circuit simulation that takes account of noise 
influences The I ff noise is understood as a stochastic process 
with a specific frequency spectrum that can be described by 
the equation 



1/f uoise sources are suitable for modeling noise influences in 
a multiplicity of teclinical and physical systems, and for sys- 
tems for estimating and forecasting events on the financial 
markets. In particular, many electronic components such as, 
for example, pn diodes and MOS field-effect transistors 35 
exhibit 1/f noise sources. 

It is possible to approximate 1/f noise sources by carrying 
out a summation of Ihe effects of many noise sources which 
each exhibit a Lorenlz- spectrum as a frequency spectrum. 
Such noise sources can be modeled, for example, by the ao 
system response of a linear time-invariant system, which can 
also be denoted as an LTI- system, at whose system input 
white noise is applied. It is disadvantageous in this mode of 
procedure for the dimension of the system of differential 
equations that are solved numerically to be excessively 45 
inflated The result of this is long computing times and a high 
storage requirement of a computer system that is used to 
simulate a system that is subject to the influence of 1/f noise. 

SUMMARY OF THE INVENTION 50 

It is accordingly an object of the invention to provide a 
method for on-demand generation of individual random num- 
bers of a sequence of random numbers of 1/f noise that over- 
comes the above-mentioned disadvantages of the prior art 53 
methods of this general type, which can be carried out quickly 
and willi a low computational outlay. It is also an object of the 
invention to specify an improved method for simulating a 
technical system that is subject to I /f noise. Finally, the aim is 
also to specify a computer system having a computer program 60 
for determining sequences of random numbers of 1/f noise 
that can he executed quickly and that lays claim to only a few 
resources of a computer system. 

In accordance with the invention, the problem of noise 
simulation in the modeling of the system to be simulated is 65 
trans formed into the problem of generating a random number 
sequence. In accordance with the invention, the correlations 
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of the random numbers are determined, and mis is used for a 
simple and accurate generation of the corresponding random 
number sequences. 
The method according to the invention for generating at 
5 least one sequence of random numbers of 1 /f noise first pro- 
vides the following steps in this case: 

a) determining a desired spectral value 3; and 

b) determining an intensity constant const, 
Characteristics of the 1/f noise to be simulated are thereby 

10 fixed. 

As a result, the number of random numbers of the 1/f noise 
that are to be generated, and a starting value for a running 
variable n used for the simulation are fixed. 

Until the desired number of elements y(n) of one or more 
15 vectors y_ of length n are calculated from 1/f-distributed ran- 
dom numbers, the invention provides the loop- type repetition 
of the following steps: 

a) increasing a current value of a running variable n by 1 ; 

b) fixing a simulation time step [t„_j; tj; 

20 c) deteraiining an element Cy of a covariance matrix C of 
dimension (nxn) according to: 



£. const -i-ty - Iff* 1 + " 'ff* 1 + I'; - '<-it* ! ~ ~ f* ' V. 
25 v 

... , n 

d) detetrnining a matrix C" 1 by inverting the covariance 
matrix C; 

e) determining a variable a in accordance with the rale 

G*$qrt( 1 /«(*.*)), 

where sqrt denotes the "square root" function, and e(n,n) 
denoting the element of the inverted covariance matrix C~ l 
indexed by (n,n); 
0 determining a (0,1 )-normally distributed random num- 
ber which forms the nth component of a vector x of 
length n; 

g) forming a variable u from the first (n-1 ) components of 
the nth row of the inverted covariance matrix C" 1 and the 
(n-1 ) elements of the vector y_ which were calculated for 
a preceding (n-1) simulation time step, specifically in 
accordance with the following rule: 




y (n _ 0 denoting the first (n-1) components of the vector 

y.» 

C~* denoting the first (n- 1 ) components of the nth row 
of the inverted covariance matrix C " 1 , and Q~* n denot- 
ing the element of the inverted covariance inalrix C~ 1 
indexed by (n,n); and 
h) calculating an element y(n) of the vector y of length n 

from l/f-distributed random numbers according to the 

following rule: 

Simulations of technical systems can be arbitrarily length- 
ened with the aid of the method according to the invention. It 
is possible for this purpose to generate additional 1/f-distrib- 
uted random numbers in a simple way when already gener- 
ated 1/f-distributed random numbers are to hand . Moreover, a 
simulation can be based on results of previously sumilated 
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time intervals. This so-called restart capability constitutes a e) forming q variables u fe in accordance with the following 

very important property for the practice of simulation. It is rule: 

precisely for I/f noise sources that this can be achieved only 

with difficulty, because random numbers that simulate a 1/f 

noise source for a specific time interval depend on random 5 y^ lxk .£> 

numbers for earlier time intervals that have already been & :» ^ 

determined numerically. The present invention also permits "* 
the use of an adaptive step size control without thereby sig- 
nificantly increasing the computation times for the simulation 

of me technical system. Such an adaptive step size control 10 y (rt _ J)Jt denotes the first (n- 1 ) components of the vector y^Uwl 
considerably increases the precision and the computation W ere calculated for a preceding simulation time step. C v j 
timeefficiency in thenumerical determination of thedynamic denotes the first (n-1) components of the nth row of the 
performance of a simulated technical system. inverted covariance matrix C" 1 andCr' denotes the element 

It is no longer necessary in the case of the method accord- - . . , nent ^„ - triY r -rj ' 



ing to the invention to prescribe the time interval to be simu- 



ofthe inverted covariance matrix CT 1 indexed by (n,u). This is 



lated. Precisely by the provision of variable step sizes, it is earned out fork I, .... q. 
also possible to perform an adaptation to current system f) calculating q elements y A ,„ which form the respective nth 
Ay na mics, and this increa ses die a ccuracy of the simulations . component ofthe vector Vj, of length n from 1 /f- distri b- 
The present invention •specifies a method for generating uted random numbers, specifically according to the fol- 
scqucuccs of I/f- distributed random numbers successively, 20 lowing rule: 
that is to say element by element. In this case, the method 
ensures rha t ca ch newly generated random number depends in y ^ m< "v* 
a correct way in the stochastic sense on the previously gen- 
erated 1 /f distributed random numbers . It is thereby possible k-i , . , . , q. 

to generate the respectively required random numbers in the 25 ^ wm ^ } of , ^ n from I/f . 

course of me numenca snnulation of a cncuit. & ^ Aq £ numte are ^ ^ iajlar 

Jne invention uses the theory of conditional probability , , . XT/ ^ rc ,~ _ ^ , , . n -J^^ n ^ 

density in order to generate a 1/f distributed random number ^antage * a NOISE matnx and specify in a snnulation the 

that correctly ensures the stochastic relationship ofthe ran- l/f noise of a system to be simulated, 

dorn number with the already generated random numbers 30 In accordance with the invention, the concept for simulat- 

requircd for preceding simulation steps. big 1/f noise is based according to the invention on the fol- 

In a particularly advantageous refmemcnt of the method lowing train of thought. The dynamic performance of a ;-:ys~ 

accordinotothemvenuoa t qscquencesofrandomnumbersof tem ^ j s expose( i to stochastic influences is modeled 

a 1/f noise are calculated simultaneously, instead of the steps adequatelv by a sl0C ha S tic process. In order to simulate such 

to be repeated m the loop-type fashion: 35 ™ 3 J . . 1 . . . . . .. 

a) determining a (0, 1 ^normally distributed random num- a s y stem S^* 1 W^ ach 13 to t M f e 
ber which forms the nth component of a vector x of numerically individual random realizations (so-called paths) 
length n; ofthe fundamental stochastic processes. Simulating systems 

b) forming a variable u from the first (n- 1 ) components of with 1 /f noise sources requires paths of stochastic integrals 0 f 
the nthrow ofthe inverted covariance matrix C 1 and the 40 the form 

(n- 1 ) elements of the vector y which were calculated for 



a preceding (n-1) simulation time step, specifically in 
accordance with the following rule: 



45 



Jo 7 



« — — 



(s)ds to be calculated numerically. Here, s is an integration 
variable and t (upper integration limit) denotes the time, 

50 

where y ( „_ 1 ) denotes the first (n- 1 ) components ofthe vector 

y f C V J denotes the first (n- 1 ) components ofthe nth row of n , (J) ^ 

IhcTinvertcd covariance matrix C"\ and C^J, denotes the 7 

clement of the inverted covariance matrix & l indexed by 

(n,n), and w 

c) calculating an element y(n) of a vectory of length n from ( s )<i s denotes a 1/f noise source, and Y(s) a stochastic process 
1/f-distributed random numbersaccording to the follow- that describes the temporal dynamics of a variable, for 
ing rule: example the electric voltage in the circuit simulation. 
r(«)-x0i)*o+ji, Using B FBJ ji$) to denote that stochastic process whose 

60 derivative (mathematically: derivative in the distribution 
the following steps are provided: sense) yie]ds ^ t J/f noise process 

d) determining a quantity q of (0,l)-normally distributed 
random numbers which form the respective last 
component of the vectors of length n, k-1 , . . . , q. It is ( ^ 
to be borne in mind in this case that the respective first 65 V" 
(n- 1 ) components ofthe vectors x* have already been 

calculated previously in die step. 
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it is possible to write the stochastic integral to be calculated as Since the individual random variables AB^^/i) represent 

increments of a Gaussian stochastic process, the random vari- 
able vector (AB FaA XO. • • • .AB/W n )> determines an n-di- 
r n mensional Gaussian-distributed random variable, and thus 

j o j (s)n^ is) ds - ^ j w> 5 ^pj^jy defined 5y its (n.dimensionai) expectation E and 

its covariance matrix C. The two variables can be calculated 
from the formulas (1 J) and (1.4) as 

Hie integra l on the right- hand side is to be interpreted as a 
Ricmann-Stieltjes integral of the stochastic process Y(s) with WW0H' s i » • » " < 3 - 5 ) 

the process B^^s) as an integrator. This integral can be 10 
approximated by a sum by decomposing the integration inter- 
val [0,t] in accordance with 05to<tj< . . . <t^t into n disjoint g : « coviAB F6M u\ ^Bn^(j)) ■ (3 6) 
subintcrvals [t„ t,_,], i=l , n: 

(1.2) 



const (-Mj-hf* % +\fj- x -r,f+ x + 
15 ]fj-r i .if' l '\t H -t i . i f 1 ),iJ^l,. 



r 



The online method according to the invention is now lo be 
specified in the form of a complete induction. 
This siun is a random variable. The dependence on the 20 Tfc e basis of the induction— and thus the starting point of 
result cd of the random experiment has been consistently the method— is the realization of a real Gaussian distribul ion 
omitted. w ith expectation 0 and variance I u =2 const ) At, \** x . 

A process B rB M whose generalized derivative has a 1/f It must be spec jfied for the purpose of inductive influence 

spectrum is known in the literature under the name of Trac- how t0 expand a rea i i2ation 0 f (AB^ 1) AB^/n-l)) 

bona! Browman Motion'. B f5A /s) is a Gaussian stochastic 25 . fl realization of £& FB jn) so as to produce overall a real- 
process and is completely characterized as such by itsexpec- nation of (AB^l), . . . .AB^n)). to order to simplify the 
ut * on notation, the already "cubed" subsequence of random num- 

£l* ra <(*)ho%e* (1.3) bers may be denoted by (y t > . . . y n _i)=:y ( „-i) r and the real- 

ization, still to be cubed, of AB^^n) may be denoted by y, r 
and by its covariance function 30 xh e problem can now be formulated as now described. 

Co^u FJ ^s)^ r ^t)^wnst(is\^ l Hti^ x -u-s^ 1 ) (1.4) Let there be given an n-dknensional zero mean Gaussian 

^ , , . „ , , random variable Zwith the covariance matrix C. Ixt the first 

The method accordmg to the invention for on-demand n _ { dements of a realization of Z be already cubed and 
generation of suitable random numbers essentially reduces m the form of a random uumber vector 

the simulahon of 1/f noise influences to the generation of . r ... ....... f _ u >u t u» „*u 

realizations of the random variables [B^t^B^^)!, , ^ 3im 18 TJ°. ^ ^ d ? tn ^ ubon *T M h * ^ 
that is to say to the generation of mcremelits of ffionai element y n mus be taken m order to complete y^to form 
Brownian Motion. a ^^on Y-<SWi>. Yn> of Z. 

The present invention permits the online generation of the A solution t0 task can be found w! f 1 <™^>» a * 
required reaUzations of the random variables: AB^ii), that 40 probability density Sfyjy^rf are considered tor y n ~with 
is to say in the course of the successive integration of the * c that y (fl _ t) already holds. Ilns variable can be 

system equations. This results in two requirements placed on calculated m the present case of a Gaussian normal distribu- 
te method: lion as 

(a) The length n of the sequence of random numbers 45 
{ABrauG), . . . ,AB FW Xn)} must remain variable dur- 
ing a simulation run. In particular, it must be possible at f(y n \ y (B _i)j 
any time to lengthen the simulation (restart capability). 
This implies the capability of the method to generate me 
additional random numbers required, for this purpose 5Q 
such that they are correlated in a correct way with the 

subsequence already generated. TTie variable Crl, is yielded inthis case from the following 

(b) Ut n be the currently attained time in the course of a notatioQ of ^ ^ rted covariailce malrix C - as: 
simulation. It must then be possible to determine the 

time interval |t„ t^, |, that is to say the step size of the 55 

next integration step, from the instantaneous system , . ^ (38 j 

dynamics— mat is to say adaptively. ^1 J rjgzli zzi 1 

Iheinveution meets both requirements by specifying a rule ~ V |£*«J 

as to how a realization of {AB FBA £\ ), . . . ♦AB^ 8JW (n)}, that is 
to say a sequence of random numbers, can be generated 50 
successively, mat is to say element by element The step size 
At,:«t -t,_j can be freely selected in this case for each new 

random number e t 1 -"*^ 

The first step is to investigate the approach for so-called £-1 e /^-n and where c' e/?. 

''conditional densities". 65 

The first step is to consider the distribution of me random 
variable vectors (AB raxj ( 1 ), . . . ,AB r ^/n)). 



j V- exp| l y(y a -n) 1 



CERTIFICATE OF CORRECTION (continued) Page 5 of 10 



US 7,308,467 B2 



8 



The variable u stands for 



v r C } 



The conditional density /(y n (y (rt _ 0 ) is therefore the prob- 
ability density of a Gaussian normal distribution with a mean 
of u and variance of 



l 



It must hold that C^,*0 in order for the above variance to 
exist. This is ensured on the basis of the following argumen- 
tation. Let C and C~ 1 have the same characteristic directions 
and inverse eigenvalues. An eigenvalue 0 of the matrix C" 1 
would therefore have an infinite variance of the random vari- 
able vector (AB^/l), . . , t AB FBJ Jn)) as a consequence. It 
may therefore be presupposed that all the eigenvalues of C" 1 
arcnonvanishing. Since the eigenvalues ofC"" 1 are non-nega- 
tive, in any case, it therefore holds that: We matrix C M is 
symmetrical and positively defined. It is possible by renam- 
ing the coordinate axes to bring this matrix from the form 
(3.8) to the following form: 



(3.10) 



This matrix is likewise symmetrical and positively defined 
per constructionem. In accordance with the Sylvester crite- 
rion for symmetrical and positively defined matrices, it fol- 
lows therefrom mat (£"'), i-£J>0, and the assertion is 
demonstrated. 

The method according to the invention reduces a simula- 
tion of 



noise sources to the generation of Gaussian-distributed ran- 
dom numbers. 

The inverted covariance matrix C~ ! (an nxn matrix) is 
required in order to generate a random number y„ that is 
correlated in the way required with an already generated 
sequence y ( „_ 0 . Strictly speaking, all that is required is 
knowledge of the nth row of this matrix, that is to say knowl- 
edge of 



v ~ *s> — an 7 



As is to be read from formula (3.6), the covariance matrix C 
depends on the decomposition of the simulation interval [0,Q 
into disjoint subintervals (step sizes) [t,_, , tj. In particular, the 
last column of C depends on t„ (doing so identically with the 
last row because of the symmetry of C), and thus on the 
current step size At^t^-t^, . 



The left-hand upper (n-1 )x(n-l) submatrix C of the nxn 
covariance matrix C is precisely the covariance matrix for a 
random number sequence of length n-1. This covariance 
(3.9) matrix had already to be determined and inverted for calcu^ 
5 lating y(n-l) or for calculating the last element (y,_,). In 
order to accelerate the method, it is therefore possible to have 
recourse to an incremental method for matrix inversion, for 
example by use of the Schur complement. 
The invention is also implemented in a method for simu- 
10 lating a technical system that is subject to 1/f noise. In this 
case, random numbers that have been determined with the aid 
of a method according to the invention are used in modeling 
and/or fixing the variables present on the input channels of the 
system. 

is Also provided is a computer system and/or a computer 
program for determining sequences of random numbers of a 
I/f noise, or for executing the other methods according to the 
invention. The invention is also implemented in a data 
medium having such a computer program, Furtbemiore, the 
20 invention is implemented in a method in which a computer 
program according to the invention is downloaded from an 
electronic data network, such as from the Internet, for 
example, on a computer connected to the data network. 
Other features which are considered as cliaracteristic for 
25 the invention are set forth in the appended claims. 

Although the invention is illustrated and described herein 
as embodied in a method for on-demand generation of indi- 
vidual random numbers of a sequence of random numbers o f 
1/f noise, it is nevertheless not intended to be limited to the 
30 details shown, since various modifications and structara) 
changes may be made therein without departing from the 
spirit of the invention and within the scope and range of 
equivalents of the claims. 
The construction and method of operation of the invention, 
35 however, together with additional objects and advantages 
thereof will be best understood from the following descrip- 
tion of specific embodiments when read in connection with 
the accompanying drawings. 

40 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a technical system to be 
simulated; 

FIG. 2 is a structogram for determining sequences of ran - 
45 dom numbers of 1/f noise; 

FIGS. 3A-3F show calculation examples for a first simu- 
lation time step; 

FIGS. 4A-4F show calculation examples for a second 
simulation time step; and 
50 FIGS. 5A-5F show calculation examples for a third simu- 
lation time step. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

55 

Referring now to the figures of the drawing in detail and 
first, particularly, to FIG. 1 thereof, there is shown a sche- 
matic diagram of a noisy system that is to be simulated. 
The system is described by a system model 1 that is 
<5o indexed as a box and describes the system behavior. The 
system behavior results from input channels 2, which are also 
designated as vector INPUT and from output channels 3, 
which are also designated as OUTPUT. Furthermore, a sys- 
tem-dictated noise is provided, which is present on uoise 
65 input channels 4 and which is also designated as a vector or as 
a matrix NOISE. The matrix NOISE is presentwhenthe noise 
is taken into account with a plurality of channels, each col- 
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uinnofthe matrix NOISE containing a vector of noise values n from 1/f-distributed random numbers. The vector y* serves 

which are present on a noise input channel. for simulating one of the noise sources for the first n simula- 

The noise on the noise input channels 4 is preferably inter- tion time steps, 

preted as a noise-dictated alteration of the system model 1. ^ch eiement of tte nt h row of NOISE is then deter- 

The behavior of the input channels 2 and of the output 5 0Q ^ ^ of equations ( 3 .7 H 3.9), from the last 

channels 3 can be desenbed by a system of differcntml equa- ^ number Qf me VBdCf x ^ ^ ., nd C) 

lions or by a system of differential-algebraic equations, so . . JLu,__ ,„ n,« /lh-_ ;„„ „.i„. 

that reliable predictions of the system behavior are possible. t0 te P™ 86 according t0 fte fo,,owuig mle - 

For each time step of the simulation o f the system shown in j^-x^c*^. 

FIG. 1 , a vector OUTPUT of the output channels 3 is calcu- 10 

Jated for a vector INPUT present on the input channels 2 and nGS 3A t0 & de P lct exemplary embodiments reproduc- 
er a vector NOISE present on the noise input channels 4. ing concrete results of calculation. 

For the simulation over a relatively long period of time, the The value of the spectral value [} is always adopted in this 

vectors INPUT, OUTPUT, NOISE are expediently specified case as 0.5. The value of the intensity const is arbitrarily 

as a matrix, a respective column k of the relevant matrix is adopted as 1.0, Three random numbers are simultaneously 

containing the values of the corresponding time series of the processed in each case in accordance with the simulation of 

relevant INPUT, OUTPUT, NOISE . no [ se sources simultaneously acting, on separate chan- 

F1G. 2 illustrates how a respective vector y* is attained nels, on the system to be simulated and which are disposed in 

which forms a column k of the matrix NOISE for the noise a vector y*, k being an integral value from 1 to 3. 

input channels 4 of the system model 1 . Each vector y* serves 20 pjQg 3A . 3 p show a calculation example for a first simu- 

for stmulatmg a noise source. ]ation ^ st [Uj M0,0.5). 

^^^^^^^ , 3A shows the cov^ance matrix C of dimension 1 xl 

sequence of computation steps is then earned out succes- & . . 0 70 *V sv, ft .u^ is t0 sav , wilh inci- 
sively for each simulation time step. re ulte u n aMl^^tion^ as 

First the current simulation time step is fixed. Equivalently res u P° na PPy bcm 1 

to this, it is also possible to fix the end of the current Simula- i .o-(-.|f l -( 1 i a5+l +ir,_ 1 -/ 1 i 0 ** l +i/ 1 -/ 1 . l i°- s+, -i/ I . 1 - 

tion time step, thereby producing the next point in time under 'i-/ * + >o^.5 M +o.5 1 - , -0=o.707iO6 . . . ; 

consideration . . . . 30 FIG . 3B shows the inverse of the covariance matrix C from 

Aftaward, the comiter nof the current simulation time step pjQ 3Af something which was performed here by uTe of a 

isincreinentedbyone. Cholesky decomposition not illustrated in more detail here. A 

The covanance matrix C of dimension (nxn) is subse- check 0 f (C C M H0.707106 . . . 0.707106 . . . ~ l ) yields the 

quently determmed accordmg to equation (3 .6). % md ^ iUuSftrate s the correctness of the value 

Ihis is followed by the step of inverting the matrix C, for for C(l I) 

example byuseofa Cholesky d^in^ition. In order tola«e pTc'sc'showsa variableaforthefct simulation stepn^l. 

theefflcency. it is also poss^le m this case to app^l to me 

inverse matrix of the previous step, for example by using J ^ 

Schur complement techniques. cr-j?rf(i ^o.707i oo . . . ), 



Next, the variable a is calculated from the formula 



40 



sqrt denoting the square root and e ( 1 , 1 ) denoting the element 
0.707106 . . . , indexed by (1,1), of the inverted covariance 

where sqrt denotes the "square root" function and where ro&fri* £ • 

e(n,n) denotes that element of the inverted covariance matrix 45 FK * 3D s^ws three values x, , x 2 , x 3 of a (0, 1 ^normally 

C ] which is indexed by (n,n). ' distributed random variable X k for one noise source to be 

" m addition, a value of a (0, 1 Hionnally distributed random simulated in each case. These values form the first element of 

variable X* is extracted and the vector x* of the normally one vector Z* of me normally distributed random num- 

distributed random numbers is thereby supplemented. The bers- The extracted random numbers have the expectation of 

extracted random number has the expectation of 0 and the so P am * me va " ance * * 

variance of 1 , This step is carried out for each noise source to FIG. 3E shows three variables for each of the three noise 

be simulated. sources to be simulated. The variable u fc is yielded in accor- 

Furthcrmore, a variable u* is formed. It is formed from the dance with formula (3.9) from the first (n- 1 ) components of 

first (n-1 ) components of the nth row of the inverted covari- the nth row of the inverted covariance matrix C* 1 as well as 

ance matrix C" 1 and from the sequence of (n-l) 1/f-distrib- 55 from the sequence of (n-I) 1/f-distributed random numbers 

uted randomTiumbers which were calculated for the preced- that were calculated for the preceding (n- 1 ) simulation time 

nig (n- 1 ) sim ulation time steps . steps . In the first simulation step, these two vectors eac h h ave 

For this purpose, the procedure is in accordance with for- hie length 0. It therefore follows for all the variables \i k in the 

mula (3 .9). This step is carried out for each noise source k to first simulation step that: u k =0. 

be simulated. 60 FIG. 3F shows three vectors y* 0 f length I of 1/f-di stri buted 

Finally, that element of the matrix NOISE is calculated random numbers that simulate the behavior of three 1 /f-dis- 

whose column index k specifies the noise source to be simu- tributed noise sources of the first simulation time step [0, 

lated and whose row index is equal to n. The current Simula- *iM0, 0.5]. The matrix NOISE is yielded from the three 

tion time step is denoted by this. The currently calculated vectors y*. The value k is an integral value of from I to 3 in this 

element r(k,n) of the matrix NOISE represents a random 63 case. Each element y k of the firs trow of NOISE is determined 

number which, together with the superjacent (n-1 ) elements on the basis of equations (3 .7}— (3.9) according to the folio w- 

of the same column k of NOISE, forms a vector & of length ing rule from the last random number x* of the associated 
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vector x and the variables ji^and a,. For example, the first 
element y,(n=l ) of the first vector y, is calculated as follows: 



10 



*,<«-l)«£t("-0* Win— 0.35 . . . 
*0.84 . . . +0.00 . . . — 0.30 . . . 

yj(n=l) and y^nH ) are calculated by analogy herewith. 

With the aid of FIGS. 4A to 4F, there is shown a calculation 
example of a second simulation time step [t, , t 2 J=[0.5, 0.75]. 
The value n for the second simulation time step is always 
equal to 2. 

FIG. 4 A shows the co variance matrix C of dimension 
(nxn)=2x2, which is required for generating one further ran- 
dom number each pernoi se source. Together with the result in 
accordance with FIG. 3F, die random number thus newly 15 
generated forms a vector y* of length 2 from 1/f-distributed 
random numbers. A vector y* is generated in this case for each 
noise source. The covariance matrix C is determined in this 
case according to equation (3.6). 

For example, this is carried out at the element C(2, 1 ) — that 
is to say with i-2 and j=l. From applying equation (3.6), 
C(2 f l)is yielded as 



12 



FIG. 4E shows three variables \x k mat were calculated from 
the inverted covariance matrix C" J in accordance with step 4B 
and from the three random numbers in accordance with step 
3F. For each noise source to be simulated, the variable \i k is 
calculated from the (n-1) first components of the nth row of 
the inverted covariance matrix C" 1 and from the sequence of 
a quantity (n- 1 ) of 1/f-distributed random numbers that were 
calculated in accordance with formula (3.9) for the preceding 
(n-1 ) simulation time steps. In the second simulation step, the 
variable U* is thus calculated from the first component of the 
second row of C 1 as well as from the first component of the 
vector y*. This is carried out, for example, with the aid of the 
value u, : 



-HI 



20 



-0.30... —1.15.. 

4.79... 

.-0.07.... 



i.o-H»i-ikl w+l +lri.i-iif fcS * l + 

h-ii.i^'-lri-i-^il"* 1 )- 
(-(0.5-0.75(° 3+1 +|0^0.75| a5+ » + 
IO.5-0.5r 5 " 1 -IO-0.51 03 * 1 )* 
-0.125 + 0.6495 ... +0- 0.3535 ... - 0.1709... 

FIG. 4B shows the inverse of the covariance matrix C from 
FIG. 4A. A check (not illustrated here) of the condition (C 
C" 1 ) yields a matrix of dimension 2x2, in the case of which 
the elements indexed by (1,1) and (2,2) are equal to 1 . The 
other elements have the value 0. 

FIG. 4C shows a variable a that is calculated from the 
inverted covariance matrix CT 1 of step 4B. The variable c is 
yielded as 

<T"sqrt(\/c(n,n)y*sqrt( 1/*(2,2))— sqrt( 1/4. 
79 . . . )-0.45 . . . ; 

sqrt denoting the square root, and e(2 , 2) denoting the element, 
indexed by (2,2) of the inverted covariance matrix C~ l from 
FIG.4B. 

FIG. 4D shows three vectors x* of independent (0, 1 )-nor- 
mally distributed random numbers, the vectors each having 
a length of 2. A (0, 1 /-normally distributed random variable x* 
is extracted per noise source to be simulated. The extracted 
random number in each case has the expectation of 0 and the 
variance of 1. The vectors x* of the normally distributed 
random numbers from FIG. 3D are thereby supplemented 
such that the vectors x* of the normally distributed random 
numbers from FIG. 4D result. 



FIG. 4F shows three vectors y* of length 2 with 1 /^distrib- 
uted random numbers that simulate the behavior of three 
1/f-distributed noise sources for the second simulation time 
30 step [tj, t 2 M0.5, o.75J. The matrix NOISE is yielded from the 
three vectors y*. The value k is an integral value of from 1 to 
3 in this case. Each element y k of the second row of NOJSE is 
determined on the basis of equations (3.7)-(3.9) according to 
the following rule from the last random number x k of the 
associated vector x and the variables u* and a. By way of 
example, the second element y t (n=2) of the first vector y x is 
calculated below: 



35 



40 



Ui(«-2hii("-2)*<J+^ l — 0.39 . . 
-0.07...-~0.10... 



» 0.45 . 



With the aid of FIGS. 5A-5F, there is shown a calculation 
example for a third simulation time step lt 2 , tjHO.75, 1 .25]. 
45 The value n is always equal to 3 during the third simulation 
time step. 

FIG. 5A shows the covariance matrix C of dimension 
(nxn)=3x3, which is required for generating one further ran- 
dom number each per noise source. Together with the result in 

50 accordance with FIG. 4F, the random number thus newly 
generated forms a vector y* of length 3 from 1/f-distributed 
random numbers. A vector y* is generated in this case for each 
noise source. The covariance matrix C is determined in this 

55 case according to equation (3,6). 

For example, this is carried out at the element C(3, 1 ) — tliat 
is to say with i=3 and j=l . From applying equation (3,6), 
C(3,l)is yielded as 



i.o.<-|f, -t 3 f™ + |/ M -f 3 |W> +1,, -Ui-i f° S4, )~ 

-{-p.5-l.2Sl"* 1 + |0-l.25| a3+, + j0.5-O.75| w+l - f0-0.75|" +l ) 
« -0 6495... +1.3975... + 0.125 - 0.6495... » 0.22... . 
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FIG. SB shows the inverse C" 1 of the co variance matrix C 
from FIG. 5A. A check (not illustrated here) of the condition 
(C C" 1 ) yields a matrix of dimension 3x3, in the case of which 
the elements indexed by (1,1), (2,2) and (3,3) are equal to 1, 
The other elements have the value 0. 5 

FIG. 5C shows a variable a that is calculated from the 
inverted covariance matrix C~ x of step SB. The variable o is 
yielded as 

<J~sqrt(]/e(rt, n)y*sqn( VeQ ,y))r-sqrt 1 0 

(1/1.75... H>. 75. . .; 

sqrt denoting me square root, and e(3 ,3) denoting the element, 
indexed by (3,3) of the inverted covariance matrix C" 1 from 
FIG. 5B. 

FIG. 5D shows three vectors x* of independent (0,1) nor- 
mally distributed random numbers, the vectors x* each having 
a length of 3 . A (0, 1 ) normally distributed random variable \ k 
is extracted per noise source to be simulated. The extracted 
random number in each case has the expectation of 0 and the 
variance of 1. The vectors x* of the normally distributed 
random numbers from FIG. 4D are thereby supplemented 
such that the vectors x* of the normally distributed random 
numbers from FIG. 5D result. 

FIG. 5E shows three variables \i k that were calculated from 2s 
the inverted covariance matrix C~ 1 in accordance with step 5b 
and from the three random numbers in accordance with step 
4/ For each noise source to be simulated, the variable a* is 
calculated from the (n-1) first components of the nth row of 
the inverted covariance matrix C" 1 and from the sequence of 30 
(n-1 ) 1/f-distributed random numbers that were calculated in 
accordance with formula (3.9) for the preceding (u-1 ) simu- 
lation time steps. In the second simulation step, the variable u* 
is thus calculated from the first two components of the third 
row of C~ l as well as from the first two components of the 35 
vector y^. This is carried out, for example, with the aid of the 
value u, : 



15 



20 



to ' 



vf c' 
...» 



-0.30... -0.31... +0.10.., 



-0.98 



1.75... 



-COO.... 



40 



45 



FIG. 5F shows three vectors y* of length 3 with 1/f-distrib- 
uted random numbers that simulate the behavior of three 
l/f-distributed noise sources for the third simulation time step 50 
[t 2 , t 3 ]=[0.75, 1 .25). The matrix NOISE is yielded from the 
three vectors y^. The value k is an integral value of from 1 to 
3 in this case. Each element y k of the third row of NOISE is 
determined on the basis of equations (3.7H39) according to 
the following rule from the last random number x^n-3) of the 55 
associated vector x and the variables \i k and o. By way of 
example, the third element y_,(n=3) of the first vector y, is 
calculated below: 



XiOrO^&Oi-*) *a+H,— -0.90 . . . • 0.75 . . . 
+0.00 , . . —0.67 . . . 

The following conditions remain to be observed for the 
purpose of concrete execution of the calculation examples 
shown. 

The numerical values shown in FIGS. 3, 4 and 5 reproduce 
intermediate and final results of the computing steps, 
described widi reference to FIG. 2, for a first, for a second and 



60 



65 



for a third simulation interval. In this case, after accurate 
numerical calculation, all values are truncated after the sec- 
ond decimal place following the point in order to be able to 
reproduce them more effectively. Consequently, when com- 
putationally reaccomplishing the exemplary embodiments, it 
is necessary to calculate further not with the intermediate 
values shown in the figures, but with the exact intermediate 
values, in order to reach the specified y vectors starting from 
the specified x vectors. 

Vectors of (0,1) normally distributed random variables are 
shown in FIGS. 3C, 4C andSC . In this case, a random variable 
represents a noise source in each case. For the sake of sim- 
plicity, it is not illustrated here how such random numbers are 
reached with the aid of the expectation of 0 and the variance 
of 1 . This is familiar to the person skilled in the art. 

We claim : 

1. A computer implemented method for numeric simula- 
tion of an electric circuit influenced by noise, comprising the 
steps of: 

numerically simulating the electric circuit using a model 
including input channels, noise input channels and out- 
put channels, the behavior of the input channels and of 
the output channels being described by a system of dif- 
ferential equations or a by a system of differential alge- 
braic equations; 

in the numerical simulation step, calculating an output 
vector for an input vector present on the input channels 
and a noise vector y of 1/f distributed random numbers 
present on the noise input channel, wherein generation 
of the noise vector y_ includes the steps of: 

determining a desired spectral value p; 

deternuning a number of the random numbers of the 1/f 
noise to be generated; 

determining an intensity constant const; 

setting a starting value for a running variable n; 

performing a loop- type repetition until a desired number of 
elements y(n)of a vector y of length n is calculated from 
1/f-distributed random numbers, by the steps of: 
increasing a current value of the running variable n by 1 ; 
setting a simulation time step [t^ ; tj; 
determining elements of a covariance matrix C of 
dimension (nxn) according to: 

V=l n 

determining an inverted covariance matrix C" 1 by inverting 

the covariance matrix C ; 
determining a variable olti accordance with 

<S-sqrt{l/e(nri) ), 

where sqrt denotes a square root function, and e(n,n) 
denotes an element of the inverted covariance matrix 
C~* indexed by (n,n); 

determining a (0, l)-nonnaily distributed random number 
which forms an nth component of a vector x of length n; 

forming a variable u from first (n- 1 ) components of an nth 
row of the inverted covariance matrix C" x and (n-1) 
elements of the vector y calculated for a preceding (n-1 ) 
simulation time step, according to: 
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where y ( „_ n denotes first (n-J) elements of the vector y, 
CrJ denotes the first (n- 1 ) components of the nth row of 
the inverted co variance matrix C~\ and C^J, denotes a 
component of the inverted co variance matrix C~ l 
indexed by (n,n); and ""5 

calculating an element y(n) of the vector y_ of length n from 
the 1/f-distributed random numbers, according to: 

outputting t he at least one sequence of random numbers of l/f 10 
noise. 

2. A method for simulating a technical system subject to l/f 
noise, which comprises the steps of: 

determining random numbers according to claim 1; and 
using the random numbers for modeling variables present u 
on input channels of the technical system. . 

3. A computer system, comprising: 

processor programmed for executing the method for deter- 
mining the sequences of random numbers of the l/f 
noise according to claim 2. 20 

4. A computer running a computer program, comprising: 
computer-executable instructions for carrying out the 

method according to claim 1 for detemuning the 
sequences of random numbers of the l/f noise. 

5. A computer running the instructions stored on a com- 25 
puter-readable data medium having the computer-executable 
instructions according to claim 4. 

6. A downloading method, which comprises the step of: 
downloading the computer program according to claim 4 

from an electronic data network onto a computer con- 
nected to the electronic data network. 

7. The method according to claim 6, which further com- 
prises using the Internet as the electronic data network. 

8. A computer system, comprising: 33 
processor programmed for executing the method for deter- 
mining the sequences of random numbers of the I/f 
noise according to claim 1. 

9. A method for simulating a technical system subject to l/f 
noise, which comprises the steps of: 40 

detennining random numbers according to claim 1; and 
using the random numbers for fixing variables present on 

input channels of the technical system; and 
simulating the technical system and outputting the result of 

the simulation. 4 5 

10. A computer implemented method for numeric simula- 
tion of an electric circuit influenced by noise, comprising the 
steps of: 

numerically simulating the electric circuit using a model 
including input channels, noise input channels and out- 50 
put channels, the behavior of the input channels and of 
the output channels being described by a system of dif- 
ferential equations or by a system of differential alge- 
braic equations; 

in the numerical simulation step, calculating an output 55 
vector for an input vector present on the input channels 
and a noise vector y of l/f distributed random numbers 
present on the noise input channel, wherein generation 
of the noise vector y includes the steps of: ^ 

detennining a desired spectral value ji; 

determining a number of the random numbers of the i/f 
noise to be generated; 

determining an intensity constant const; 

setting a starting value for a running variable n; $5 

calculating q sequences of the random numbers of the l/f 
noise simultaneously, by performing loop-type repeti- 
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tions until a desired number of elements y kw „ of a vector 
y of length n is calculated from l/f-distributed random 
numbers, by the steps of: 

increasing a current value of the running variable n by 1 ; 
setting a simulation time step [t„_, ; tj; 
determining elements of a co variance matrix C of 
dimension (nxn) according to: 

IJ m \ 9 . . ,n 

determining an inverted covariance matrix C" 1 by inverting 

the covariance matrix C; 
determining a variable a in accordance with 

a»sqrtl\fe(nn)), 

where sqrt denotes a square root function, and e(n,n) 
denotes an element of the inverted covariance matrix 
C~ l indexed by (n,n); 

detennining a quantity q of (0, 1 )-normally distributed ran- 
dom numbers x k n which form a respective last compo- 
nent of vectors x* of length n, where k=l, . . . , q, 

forming q variables u fc according to: 



^T— 

where y (n _ t ^ k denotes first (n- 1 ) elements 0 f the vector y^ 
that were calculated for a preceding simulation time 
step, Cr,| denotes the first (n-I) 
components of the nth row of the inverted covanance 
matrix C" 1 , and C^J, denotes the elem ent 0 f the inverted 
covariance matnx CT l indexed by (n,n), where 
k=l, . . . , q; and 
calculating q elements y* n which form a respective nth 
component of the vector of length n from 1 /f- distrib- 
uted random numbers, according to: 

where k-1, . . . , q; 

outputting at least one of the q sequences of random num- 
bers of l/f noise. 

11. A method for simulating a technical system subject to 
l/f noise, which comprises the steps of: 

determining random numbers according to claim 10; and 
using the random numbers for modeling variables present 
on input channels of the technical system. 

12. A computer system, comprising: 
processor programmed for executing the method for deter- 
mining the sequences of random numbers of the l/f 
noise according to claim 11. 

13. A computer running a computer program, comprising: 
computer-executable instructions for carrying out the 

method according to claim 10 for detennining the 
sequences of random numbers of the l/f noise. 

14. A computer running the instructions stored on a com- 
puter-readable data medium having the computer-executable 
instructions according to claim 13. 

15. A downloading method, which comprises the step of: 
downloading the computer program according to claim 13 

from an electronic data network onto a computer con- 
nected to the electronic data network. 

16. The method according to claim IS, which further com- 
prises using the Internet as the electronic data network. 
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17. A computer system, comprising: 

processor programmed for executing the method for deter- 
mining the sequences of random numbers of the 1/f 
noise according to claim 10. 

18. A method for simulating a technical system subject to 
1/f noise, which comprises the steps of: 
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determining random numbers according to claim 10; and 
using the random numbers for fixing variables present on 

input channels of the technical system; and 
simulating the technical system and outputting the result of 

the simulation. 

* * # * * 



